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Abstract 

Based on symmetry consideration of migration and shape deformations, we formulate phenomenologically the dynamics 
of cell crawling in two dimensions. Forces are introduced to change the cell shape. The shape deformations induce 
migration of the cell on a substrate. For time-independent forces we show that not only a stationary motion but also a 
limit cycle oscillation of the migration velocity and the shape occurs as a result of nonlinear coupling between different 
deformation modes. Time-dependent forces are generated in a stochastic manner by utilizing the so-called coherence 
resonance of an excitable system. The present coarse-grained model has a flexibility that it can be applied, e.g., both to 
keratocyte cells and to Dictyostelium cells, which exhibit quite different dynamics from each other. The key factors for 
the motile behavior inherent in each cell type are identified in our model. 
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1. Introduction 

Eukaryotic cell crawling has attracted much attention 
recently from the view point of nonlinear science and non¬ 
equilibrium statistical physics. One of the characteristic 
features is that the symmetry is spontaneously broken to 
cause the front-rear asymmetry when the cell migrates in 
contrast to bacterias which swim by rotary motion of flag¬ 
ella and hence are inherently asymmetric. The dynamics 
of eukaryotic cells involves the mechanical forces between 
cell membrane and substrate, and biochemical reaction of 
active molecules inside a cell. 

Study of cell crawling on substrates has began rather 
recently compared to that of swimming bacterias. The 
latter has a long history of hydrodynamical approach in 
the limit of low Reynolds number [ll-y. Shape deforma¬ 
tion of crawling keratocyte cells has been analyzed exper¬ 
imentally [^. Classification of morphology of motile cells, 
correlations between shape deformation and migration of 
Dictyostelium cells and other living cells have also been 
investigated [^-0. Recent advanced experimental tech¬ 
niques have enabled us to measure the spatial distribu¬ 
tion of traction forces exerted by a migrating cell on sub¬ 
strates [IMl and the concentration distribution of active 
molecules which involve cell motility [l3| . 

Plasma membrane protrusion caused by actin polymer¬ 
ization in the cell interior is the essential mechanism of cell 
crawling. In the early 1990’s, DiMilla et al investigated 
persistent migration of tissue cells such as fibroblasts by a 
mathematical model which is essentially one-dimensional 


and incorporates cytoskeletal force generation, cell polar¬ 
ization, and dynamic adhesion [l5| . Theoretical studies 
of cell crawling taking into account shape deformations 
have started recently. Modeling of cell crawling employ¬ 
ing reaction-diffusion mechanism inside a cell or on a cell 
boundary and the interaction between the chemical com¬ 
ponents and the cell membrane has been proposed in two 
dimensions 16 -1^. A phase field model for cell shape 


coupled with the polarization field of actin filaments has 
also been proposed for the cell motility. The crawling 
dynamics of keratocyte cells including oscillatory straight 
motion and bipedal motion [l^ has been investigated by 
taking account of the density of adhesion bonds and trac¬ 
tion forces 2^, 211. A similar model in terms of the phase 
field has been studied in which the reaction-diffusion dy¬ 
namics inside a Dictyostelium cell is assumed to be ex¬ 
citable and therefore this model is capable of investigating 
non-stationary motion 1^. Modeling of amoeboidal cell 


crawlin g h as also been formulated by an oscillatory dy¬ 
namics where irregularity appears as spatio-temporal 
chaos [23|. These are models in two dimensions. In a 
slightly different approach, theory of active gels has been 
applied to stationary amoeba motion in one dimension to 
make a connection between the migration velocity and the 
distribution of active stress or myosin molecules 24 l25|. 


Motility of active droplets in which active stress is gener¬ 
ated has been studied numerically both in two and three 
dimensions to show that a bifurcation from a motionless 
state to a migrating state with shape deformations occurs 
due to spontaneous symmetry breakin g of the polarity in¬ 
version in the absence of treadmilling [26| |. 
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It is mentioned briefly that there is another mechanism 
of cell motility due to plasma membrane blebbing 271. 
This is not restricted to migration on substrate. Bleb¬ 
bing is initiated by local disruption of membrane-action 
cortex and internal hydrostatic pressure. It is of impor¬ 
tance to note that actin polymerization is not involved 
in the initial bleb expansion. Cell motility by blebbing 
in three-dimensional environments has been investigated 
theoretically 281. 

In the present paper, we study cell crawling under a 
homogeneous environment based on a phenomenological 
model in terms of migration velocity and shape defor¬ 
mations. A set of time-evolution equations is derived 
based on symmetry consideration in the same spirit as the 
derivation of equations for deformable self-propelled par¬ 
ticles [Hilo . To make shape deformations, we introduce 
forces which act on the cell perimeter. We consider time- 
independent and time-dependent forces separately. The 
case of constant forces is regarded as a model of coher¬ 
ent motions, e.g., of keratocyte cells whereas the time- 
dependent forces are applied to motility of Dictyostelium 
cells. In the experiments of Dictyostelium cells, morpho¬ 
logical change occurs repeatedly but it is not precisely pe¬ 
riodic. To realize this behavior, we utilize the so-called 
coherence resonance which generates spike excitation of 
chemical components repeatedly in an excitable system 
when noise is added appropriately [3ll |. 

All the previous models mentioned above are con¬ 
structed to apply to a specific system such as fish kera¬ 
tocyte cells and Dictyostelium cells. A steadily migrating 
keratocyte cell is elongated perpendicularly to the veloc¬ 
ity direction whereas a Dictyostelium cell in a starved 
condition has a tendency to elongate parallel to the migra¬ 
tion direction Q. Our model, though simple, has an ad¬ 
vantageous feature that it is applicable to non-stationary 
motion of a crawling cell with general shape deformations 
by choosing appropriately the parameters. 

One of the basic assumptions of our model in terms of 
the center of mass and the cell boundary is that all other 
degrees freedom involving migration relax rapidly. If this 
is not the case, we need to add other relevant degrees of 
freedom as dynamical variables. 

In the next section (section 2) we start with descrip¬ 
tion of our model system. The case of constant forces is 
analyzed in section 3 where we obtain limit cycle oscilla¬ 
tions of migration velocity and shape deformations. This 
is compared with the phase field model of keratocyte cells. 
Numerical results for the time-dependent forces are shown 
in section 4 and are compared qualitatively with the mo¬ 
tions of Dictyostelium cells. Discussion is given in section 
5. 


2. Model for cell crawling 


We introduce the model of cell crawling in two dimen¬ 
sions in terms of the velocity of the center of mass Vi, and 


the deformation tensors 

'^SijUijkj (1) 

S' 

-K2Sij + boiviVj - ^VkVk) + Fj;^\t)[2) 

-KsUijk + do[viVjVk 

+ 6jkVr + dkiVj)] + T;^/2(t),(3) 

where the repeated indices imply summation. 

The tensors Sij and Uijk are defined as follows. Defor¬ 
mations of a cell around a circular shape with radius i?o 
are written as 


Vk = 
dSij 

dt 

dUijk _ 

dt 


R{(j),t) = i?o(l -I- 5R{(j),t)) , 
where 

OO 

6R{<P,t)= Y. 

n— — oo 


( 4 ) 

( 5 ) 


Since uniform expansion and contraction of a circular cell 
are prohibited and the translational motion of the cell has 
been incorporated in the variable Vk, the modes cq and c±i 
should be removed from the Fourier series ([S]). The defor¬ 
mation tensors are given in terms of the Fourier coefficients 


by M 

Sll = C2 -I- C_2 = 202 COS202, (6) 

<S'i2 = S'21 = i(c2 — c_2) = 202 sin 202 , ( 7 ) 

S 22 = —Sll, (8) 

C/111 = —C/122 = —C/212 = —C/221 = kF+, ( 9 ) 

C/222 = —C/112 = —C/121 = —C/211 = —W-, ( 10 ) 

where 

W+ = C 3 + C _3 = 203 cos 36»3, (11) 

W- = i{c 3 - C- 3 ) = 2o 3 sin36»3, (12) 


with positive 02 and 03. 

The coefficients K 2 and K 3 are positive while the sign 
of 7 in eq. 0 will be fixed later. Here, for simplicity, 
we ignore other nonlinear couplings such as UijkVk and 

Sii Vk SjkUi SkiVj {Vji/Skn SjkSin SkiSjn} 

[2^ but consider the coupling only with the velocity as 
eqs. @ and ([21) since those terms are expected to be 
mostly relevant to the correlation between the elongation 
direction and the migration direction as shown below. 

Equation 0 implies that there is no inertia term and 
that the cell does not migrate if it is circular since we con¬ 
sider a deformation-induced migration. In our previous 
studies of migration-induced deformations H 351, equa¬ 
tion of motion of the center of mass was derived, which 
takes the following form 

*?(^z) ‘^k T ClSkjVj -\- 7 SijUijk, (48) 
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where ki, g{> 0), a and 7 ' are constants. Equation ([T]) is a 
special case of eq. (IT^ . In fact, when m is negative, that 
is, migration is passive, one may ignore the g-term and in 
the limit of |ki| » 1 , the solution of eq. (USD is given by 

Vk — 7 [^1 “f ~ SijUijk- (14) 

Ki 

It is also noted here that the deformation-induced migra¬ 
tion in the form of eq. o has been introduced in a dif¬ 
ferent context of cell motility [s^. In our model, the coef¬ 
ficients K 2 and K 3 are positive assuming that the circular 
shape is always stable when the deformation forces are ab¬ 
sent. The front-rear asymmetry is produced when the cell 
starts to migrate. 

Deformations are caused by the internally created forces 
given by the terms and F-p in eqs. ([ 2 ]) and Q, 
respectively, which are called deformation forces. In the 
present model, these are force moments defined on the cell 
boundary. Let us suppose a force acting on the 

position a on the boundary in a two-dimensional cell. The 
moments are defined by 

Fp\t) = J daxifj{a,t), (15) 

F^pit) = J daXiXjfkia,t), (16) 

where the integral runs over the cell boundary. In experi¬ 
ments of crawling cells, the distribution of traction forces 
between the cell membrane and the substrate has been 
measured [l^ . To make a link between the traction forces 
and the deformation force tensors, one needs further in¬ 
formation inside the cell such as the viscoelastic proper¬ 
ties. In the present simple model approach, we do not 

make such a systematic reduction of degrees of freedom 

but simply assume that these force tensors have the same 
symmetry as that of the deformation tensors given by eqs. 
® - (HU)- For example, we redefine F^p — iSij/2)Fp'’ as 
F }-' to make the tensor traceless. In a manner consistent 

•‘J _ _ _ 

with the expressions ®, ®, (dH) and (HU), we write the 
force tensors as 

Etf = 5^2)(i)cos20g\ (17) 

F^P = cos (18) 


The other is to equate it with the velocity direction, i.e.. 


9 ( 2 ) 

^11 

9 ( 2 ) 

^12 

9 ( 3 ) 

'-'ill 

Cl(2) 



0. 



(23) 

(24) 

(25) 

(26) 


where the angle 9 is related with the velocity of the center 
of mass as 


vi = vcos9, (27) 

V2 = usin 0 , (28) 

with the amplitude 1; > 0 . 

In the present paper, we mainly consider the case given 
by eqs. (fro]) - (l^ . In this case, the deformation forces are 


given by 



^(2) 

^11 

= (5^^^^)+ 6 ( 0 ) COS 2(02-f 772(0), 

( 29 ) 

^(2) 

^12 

= ( 5 ^^ 0 O + 6 ( 0 ) sin 2(02+ 772(0). 

( 30 ) 

pi3) 

^111 

= (S^^OO+ 6 ( 0 ) cos 3(03 -b 773(f)), 

( 31 ) 

^( 3 ) 

^112 

= (fl^^OO + 6 ( 0 ) sin 3(03 -b 773(0), 

( 32 ) 


where we have added random noises 72 , Vs, ^2 and ^3 both 
to the amplitudes and to the phases. We assume that these 
noises obey the Gaussian statistics with zero means and 
are delta-correlated in time; 


< ?7n > 

< > 


<^n>=0 (n = 2,3) (33) 

al6nmS{t-t') (n = 2,3) (34) 
CSnmS{t-t') (n = 2,3) (35) 
0 (n = 2,3) (36) 


with (T„ and the magnitudes of the noises. We do not 
add a random force to eq. o but explore the consequence 
that randomness of a migrating cell is generated by the 
stochasticity of the shape dynamics. 

Equation o can be written as 

■^1 = 7 [('S'ii ~ <522)^111 — 2S12U222] 

= 7[2s2S3Cos(303-202)], (37) 

V2 = 7[('S'22 — Sii)U222 — 2S'i2Giii] 

= 7[2s2S3sin(303-202)], (38) 


where the amplitudes g^^\t) and ^^^^(t) will be specified where S 2 — 2a2 > 0 and S 3 — 203 > 0. Using eqs. ([27]) 
later. There are basically two choices of the phase ([211)) we obtain 


One is to equate it with the deformation directions 
tently with eqs. ®, ®, (fTTl) and (IT^: 

consis- 

V = 2|7|S2S3, 

0 = 303 - 202 -4'„, 

(39) 

(40) 

= 02, 

(19) 

where 


eg2^ = 6-f, 

( 20 ) 

= 0 , 

(41) 

©S = 6 , 

( 21 ) 

for 7 > 0 and 


©^21 = 6 + f. 

( 22 ) 

'^v = tt. 

(42) 
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for 7 < 0 , 

Equations © and ([3]) are written in terms of the ampli¬ 
tudes and angles of deformations as 

^ = -K 2 S 2 -h ^^cos(6d23 + 2 ^'„) 

+ (43) 

§ = + 2 *,)+ (44) 

at 4s2 S 2 

^ = -K3S3 + cos(6023 + 

+ 9^^Ht)+^3{t), ( 45 ) 

^ = -^^sin(6023 + 3«'„) -h — —r]3{t), (46) 

at 12 s 3 S 3 

where the relation (Hill) has been used. Since the system is 
isotropic, the relative angle defined by 

023 = 02 — 03, (47) 

enters in the time-evolution equations. We have ignored 
the terms nonlinear in the noises employing the approx¬ 
imations such that (cos 2772 ) —t 1, ( 005373 ) —?> 1 and 
^3(t)(sin373) -)> 0. 

Now we discuss the amplitudes of the deformation forces 
9 ^‘^\t) and g^^^t). When our model is applied to ker- 
atocyte cells which migrate coherently without irregu¬ 
lar shape deformations, the deformation forces are as¬ 
sumed to be constant. On the other hand, experiments of 
Dictyostelium cells show that the shape changes repeat¬ 
edly such that a cell elon gate s first and then the third 
mode deformation follows |13l |. Therefore, the deforma¬ 
tion forces in this case are time-dependent with a suitable 
phase difference between the modes. This means that if 
the forces are approximated by periodic functions with pe¬ 
riod T such that 

= g^^\t + T), (48) 

g(3)(i_$) = + T -4>), (49) 


the phase difference $ should be 0 < $ < tt. 

Another possibility is that generation of the deformation 
forces has an excitable origin. As mentioned in Introduc¬ 
tion, Sawai and his coworkers have found experimentally 
that motility of a Dictyostelium cell is coupled with chem¬ 
ical reactions inside the cell and they have analyzed the 
dynamics by a set of excitable reaction diffusion equations 
with noises 1^ . In the present model approach, we utilize 
the so-called coherence resonance 3l| to creat stochastic 
but almost periodic deformation forces. The phase differ¬ 
ence between the modes as eq. (091) is also considered as 
shown in section 4. 

From the solutions of eqs. dSH), dMl), (031), (gH), (051)- 
and ( 00 ). the motion of a migrating cell is represented as 


xit) = Xcm{t) + r{t)cos{a), (50) 

y{t) = Vcniit) +r{t)sm{a), (51) 


(a) 


(b) 
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Figure 1: Stationary shapes of a crawling cell subjected to constant 
deformation forces for (a) Rq = 0.5, S 2 = 0.1, S 3 = 0.3 and 62 = 
63 = 0, and (b) Rq = 0.5, S 2 = 0.3, S 3 = 0.1, 02 = 0 and O 3 = —tt/O. 
When 7 < 0, the cell migrates to the left in (a) and upward in (b). 


where 0 < a < 27r and 

r{t) = i?o (^1 -I- S 2 cos 2(a - 02 ) + S 3 cos 3(a - 6 * 3 ) j, (52) 

with (xcm,ycm) the position of the center of mass of the 
cell. 


3. Numerical results for constant deformation 
forces 


Hereafter throughout present paper, we consider the 
case 7 < 0 and = tt (eq. HU) since this condition 
is found to correspond to motility of keratocyte cells and 
Dictyostelium cells. 

In this section, we analyze the case that the deformation 
forces are independent of time and noises are absent. From 
eqs. (03|). (031), (03) and (H^ ignoring the noise terms, we 
obtain 

^ = -K2S2 + cos(66»23) + (53) 

^ = -K 3 S 3 - cos(66»23) + (54) 

at 4 


d023 

dt 

where 


-C'sin(66»23), 


bov^ ^ dov^ 
4^ 1 ^' 


(55) 


(56) 


We require that the amplitudes S 2 and S 3 are positive. 
This condition is fulfilled for the stationary solutions of 
eqs. dMl and (IMl) with eq. (15^ if I 7 I is not very large. 
We put 02 = 0 in eq. (15^ without loss of generality. That 
is, the long axis of an elliptical deformation is parallel to 
the X axis. If C is positive which is realized for bg > 0 and 
do > 0 , 023 = 0 or 7 r /3 are a stable solution. The former 
gives us 03 = 0 and 0 = —tt. The stationary shape in this 
case is displayed in Fig. (Ha). The cell is migrating to the 
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left at a constant speed. The other solution 023 = tt/S 
gives us 03 = —tt/S and 0 = —2tt. Therefore the motion 
and the shape in this case is just mirror symmetric with 
respect to the y axis of Fig. [DJa). When C is negative as 
is realized by the condition that 6 o < 0 and (Iq < 0 , the 
stable solution of eq. (1551) is given by 023 = ±7r/6. The 
shape for 03 = —tt/G is shown in Fig. [IJb). The angle 
of the velocity direction is given by 0 = — 37r/2 indicating 
that the cell migrates upward. The other case 03 = tt/G 
is just mirror-symmetric with respect to the x axis. The 
steady motion of keratocyte cells is qualitatively expressed 
as in Fig. HTb) whereas a shape of Dictyostelium cells more 
resembles to that in Fig. [DJa). When 7 is positive, the 
migration direction is reversed without any influence for 
the cell shape. 

We investigate the dynamics of a cell in Fig. [IJb) for 
negative values of bo and do in further detail. Let us in¬ 
troduce 823 through the relation that 

^23 = w + <523- (57) 

G 

Equations (1551) . (I5il) and (1551) hold after replacement of 
023 by ( 523 , bo by | 6 o| and do by |do|- Note that 823 = 0 is a 
stable solution. The nonlinear terms in eqs. (1531) and (1551) 
have an interesting structure. Note that v is proportional 
both to S 2 and to S 3 . If S 3 is large, it makes S 2 increase 
in eq. (insD. This, in turn, tends to decrease S 3 by the sec¬ 
ond term in eq. (1541) . However, when the time-variation 
of S 3 is sufficiently slow, i.e., K 3 <C /C 2 , there is a delay of 
this evolution and an overshooting can occur, which make 
an oscillatory dynamics of S 2 and S 3 for constant values 
of and We have verified numerically that this 

is indeed the case. A supercritical bifurcation occurs for 
7 = —1.0, K 2 = 1, = do = —3.5 and g^^^ = g^^'^ = 0.1 at 

around ~ 0.125 below which the stationary solution of 
eqs. dSSI) and (1551) looses stability and a limit cycle oscil¬ 
lation emerges. Figure [2] displays snapshots of a migrating 
cell for K 3 = 0.12 in the interval 175 < t < 200. The corre¬ 
sponding time-variations of v, S 2 and S 3 are shown in Fig. 
[3] The period of oscillation T is about T = 17. 

The oscillatory straight motion in Fig. [3]should be com¬ 
pared with the one obtained in the phase field model of 
keratocyte cells Figure 2 in ref. indicates that 
the aspect ratio d of a cell takes a maximum just after the 
velocity becomes a maximum. The asymmetry parameter 
C exhibits approximately an anti-phase oscillation with the 
aspect ratio. See ref. [20j for the definitions of h and (^. 
If we assume the correspondence d — 1 o S 2 and C O S 3 , 
which is valid for 82,83 <C 1 , the dynamical behavior in 
Fig. [3] is quite similar to the one in ref. [ 2 ^. Only the 
difference is that S 3 > S 2 in our model whereas d — 1 > ^ 
in the phase field model so that the shape of an oscillating 
cell looks slightly different. It is also noted that the phase 
difference between S 2 and S 3 depends on the model pa¬ 
rameters. For example, for 7 = —1.0, K 2 = 1.0, K 3 = 0.05, 
bo = —4.0, do = —3.0, = 0.12 and g^^'> = O.OG, the am¬ 

plitude 82 reaches the maximum before S 3 takes its mini- 



Figure 2: Snapshots of a cell undergoing an oscillatory rectilinear 
motion from the top (at t = 175) to the bottom (at t = 200). The 
parameters are set as 7 = —1.0, k .2 = 1, K 3 = 0.12, 60 = do = —3.5 
and = 0.1. The symbols x indicate the trajectory of the 

center of mass. 



Figure 3: Time-evolution of v (red and solid line), S 2 (green and 
broken line) and S 3 (blue and dotted line). The parameters are the 
same as those in Fig. [2] 


mum value. 

The oscillatory rectilinear motion occurs also for the pa¬ 
rameters of Fig. [Ha) where the cell elongates parallel to 
the migration direction. It is readily shown that when 
7 < 0 , and bo and do are positive, the same argument of 
oscillation as the above can be applied to eqs. (1551) and 
dMD with 023 = ^ 23 . Therefore the phase difference be¬ 
tween S 2 and S 3 is approximately equal to tt. However, 
experiments of Dictyostelium cells indicate that a migra¬ 
tion cell tends to first elongate along the velocity direction 
and, shortly after this, grows a Y-shape deformation 13|. 
This fact is not in accord with the phase difference between 
S 2 and S 3 mentioned above. To achieve the experimentally 
observed sequence of the shape change, we need to intro¬ 
duce time-dependent deformation forces as described in 
the next section. 

Finally we make a remark that if we add the fourth 
mode, c±i{t) in the time-evolution equations, a circular 
motion appears as a stable solution after destabilization 
of a straight motion, but, other oscillatory motions such 
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as a bipedal motion and a zig-zag motion have not been 
obtained at present. 


0.25 


4. Numerical simulations for time-dependent de¬ 
formation forces 


In this section, we put 7 < 0, 69 > 0 and do > 0 to in¬ 
vestigate the dynamics of a migrating cell as shown in Fig. 
[TJa) under the influence of the time-dependent deforma¬ 
tion forces and noises. First, we need to create the forces 
and Experiments of Dictyostelium cells 

indicate that the force dipole and the force quadrupole 
are time-dependent but not precisely periodic [l^. The 
phase difference between these two is another important 
parameter in constructing a model system. Furthermore, 
Sawai and his coworkers have interpreted the motility of 
Dictyostelium cells as an excitable system 1^. In order to 
take this into account and to realize an appropriate phase 
difference in a stochastic excitable system, we employ the 
coherence resonance 31|. That is, we consider the follow¬ 
ing excitable system; 


dw 


T 


dt 

dz 


dt 


a{w - -w^ - z), 
a{w -I- a -I- Cz), 


(58) 

(59) 


where tr is a positive constant which is needed to adjust 
the average period of the deformation forces. The noise Cz 
is Gaussian with 


< €z > = 0 , (60) 

< ez(.t)ezit') > = d^5{t-t'). (61) 

where d is the noise intensity. 

Here, we briefly summarize our method of numerical 
computation of the stochastic equations eqs. (liHl) . (HI) . 
(03, (03 and (03, all of which take the following form 

^ = f{X) + h{X)z, (62) 

where the Gaussian noise 2 satisfies 


< z > = 0, (63) 

< z{t)z{t') > = Dx6(t — t'). (64) 

with the noise intensity Dx- We employ the standard 
Euler-Maruyama scheme [s^ for the discretized version as 

X„+i - = At/(X„) + {Aty^^DxGnhiXn), (65) 

where is Gaussian random number with < Gn >= 0 
and < >= 1, and At the time increment. For sim¬ 

plicity, we replace G„ by another random number '/SG!^ 
where G'^ is distributed uniformly between -1 and 1. It 
has been proved mathematically that numeral accuracy is 
unaltered by this replacement for sufficiently small At [33| . 

In eq. (1^ . we put d = 1/(10-^^). When a > 1, the 
solution w = —a, and z = —a + a^/3 is linearly stable 




Figure 5: Time-evolution of gG) (red full line) and (green broken 
line). 


and the system is excitable when r is very small. We have 
examined the dynamics in the regions 1.01 < a < 1.1 and 
10 < 1/r < 100. For larger values of a and for smaller 
values of l/r, the time interval between two adjacent exci¬ 
tations increases, but the duration time of each excitation 
is insensitive to these parameters. 

The deformation forces are assumed to be given by 

= g^^\w/3 - 2V2zli + 5 ^^ ( 66 ) 

g(3) = gf\2V2w/?. + z/3 + g^i^), (67) 

where we put g^^^ = 0.07, g^^ = 2, = 0.05 and g^^ = 

2.5 unless stated otherwise. The behavior of and 
for r = 0.01, a = 1.05 and cr = 0.3 is shown in Figs. 0] 
and [5l Note that g^^'> increases first while g^^'> remains 
small and almost constant, and then both exhibit a rapid 
increase. After reaching the maximum, the forces decrease 
rapidly and the process is repeated. Note also that this 
dynamics is noise-controlled and is not strictly periodic as 
can be seen in Fig. [SJ 

We have solved numerically eqs. ( 07 ]), ( 03 ), ( 03 , ( 03 , 
( 03 , and (03 with the forces (03 and 03 . As men¬ 
tioned above, the Euler-Maruyama scheme with the time 
increment At = 10“^ is used in the presence of noises. Fig¬ 
ure 13 displays a migrating cell without noises Un = ?)« = 0 
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Figure 6: Snapshots of a cell migrating to the left in the absence of 
noises at t = 966.25, t = 968.75, t = 971.25, t = 973.75, t = 975.25 
and t = 977.75. The straight line is the trajectory of the center of 
mass. 


for the parameters 7 = —4, Rq = 0.1, bo = do = 1.0, 
K 2 = K 3 = 1-0. An almost circular cell first elongates 
in the direction of the migration velocity, then the third 
mode of deformation grows so that the cell is accelerated. 
As the deformations become small, the cell gets slowdown 
and relaxes to the circular shape. 

Numerical simulations in the presence of random noises 
have been carried out both for the linear case bg = do — 0 
and for the nonlinear case bo = do = 2. Other parameters 
are fixed as 7 = — 10 , K 2 = K 3 = 1 . 0 , €amp = 'J^On =■ 0.02 
and €ang = = 0.2 (n = 2,3). Figure [7] shows 

the results for the linear case. The snapshots of the mi¬ 
grating cell with Ro = 1.0 are shown in Fig. [UJa). The 
time-evolutions of the angles 0 , O 2 and O 3 are depicted in 
Fig-EDb) whereas those of S 2 , S 3 , v and 6 are in Fig. [TDc). 
The discretizing version of 0 = 303 — 202 as eq. (1551) di¬ 
vided by At is used in the plot of 0. The corresponding 
hgures for the nonlinear case are given in Figs. ISa), (b) 
and (c), respectively. The radius of the cell in Fig. |SDa) is 
chosen as Ro = 2.0 because the liner size of this figure is 
about twice as large as that of Fig. EDa). 

Now we compare the two cases in detail. There are dis¬ 
tinct differences in the dynamics between the linear and 
nonlinear cases. The trajectory of the linear case in Fig. 
Efa) is much more random compared with that in the non¬ 
linear case in Fig. |5Da). In other words, a cell having the 
nonlinearity exhibits orientational persistence in its migra¬ 
tion. In the present set of parameters, the nonlinear terms 
make the migration velocity parallel to the elongation di¬ 
rection of the cell whereas the deformation directions and 
the migration direction are determined solely by the ran¬ 
dom noises as in eqs. (HI) and (H5)) when the nonlinear 
terms are absent. One unexpected behavior is that a cell 
reverses sometimes the direction of migration in the non¬ 
linear case as can be seen around at a; = 38 and y = —13 
in Fig. ISa) although this does not occur for weaker noise 
intensities. 


One may note that the time change of the migration 
direction 0 in Fig. EDb) looks similar to that in Fig. [S^b) 
apart from their amplitude. This is because these data 
were obtained by the same sequences of random noises. 
We have verified the random motion in the linear case and 
the persistent motion in the nonlinear case by other several 
independent runs of numerical simulations. 

Although the motion in Fig. EDa) looks like a Brownian 
motion, the center-of-mass motion and the shape change 
are not totally independent of each other because of the 
relations dSSl) and (pni) . It is evident in Fig. El^b) that, 
in the linear case, 0 exhibits a time-variation similar to 
03 but there is no apparent correlation between 0 and 02 . 
In contrast, these three angles change in an almost syn¬ 
chronized manner in the nonlinear case as can be seen in 
Fig. [HDb). Note that the variation of 0 is in the range of 
0 < 0 < 27r (modulus 27r) so that the trajectory is curved 
substantially in the linear case whereas it is smaller in the 
nonlinear case as 0 < 0 < tt (modulus 27r) so that the 
trajectory does not exhibits much wandering. 

The migration velocity is determined by the amplitudes 
S 2 and S 3 as eq. (ESI). Note that the magnitude of the 
migration velocity is not much different between the lin¬ 
ear and nonlinear cases. Therefore, the longer persistence 
length in the nonlinear case is not due to velocity differ¬ 
ence. The two deformations S 2 and S 3 take a maximum 
almost simultaneously in the linear case as in Fig. EJc) 
because of the time-dependence of the forces as shown in 
Fig. El In the nonlinear case depicted in Fig. EKc), the 
amplitude S 2 becomes maximum with a slightly delay com¬ 
pared with S3. The value of S3 is much smaller than that 
of S 2 - This is a consequence of the fact that the nonlinear 
term in eq. (l43l) has a plus sign but that in eq. (l45l) has 
a minus sign for 7 < 0. The magnitude of fluctuations of 
0 indicates an interesting time-variation as shown in Figs. 
E{c) and[5Kc). In the linear case, it becomes large just 
at the instance that the velocity (and S2 and S3) starts 
increasing. In the time region when these quantities are 
decreasing, the fluctuations of 0 are also small. The main 
reason of this is due to the multiplicative noises for the 
angles in eqs. (PI) and (PI). Because of a delay in the 
evolutions of the deformations S 2 and S 3 relative to the 
deformation forces, the cell is not deformed substantially 
even when the forces are large and vice verse in the linear 
case. In the nonlinear case, 0 fluctuates strongly during 
the growth of the velocity whereas appreciable decrease 
in fluctuations is not seen in the decreasing of the veloc¬ 
ity since the nonlinearity masks the multiplicative effect of 
noises. 

We have assumed in this section that the origin of the 
deformation forces is excitability. The case of oscillating 
deformation forces will be described shortly in the next 
section. 

Before closing this section, we compare our results with 
experiments of Dictyostelium cells in a vegetative con¬ 
dition [ 1 ^ to check the consistency of the characteristic 
space-time scale. The characteristic period of the force 
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Figure 7: (a) Trajectory of the center of mass (green line) and the 
shape change (red lines) in a cell migration for the linear model. 
The snapshots are taken for 812.5 < t < 987.5 with the time-interval 
St = 25. The cell starts to migrate around at a: = —4.0 and y = 11 at 
t = 800. First it moves downward and then upward to a: = —11.4 and 
y = 14.9 at t = 1000. (b) Time-evolution ol A = 9/{27r) (red line), 
B = 02/(27r) (green line), and C = 9^/{2n) (blue line) for the linear 
model. It is evident that the velocity direction 9 exhibits a correlated 
time-change with the angle of the third mode of deformations 9$. 
Such a correlation is less clear in the angle 92 of the second mode, 
(c) Time-evolution of v (red line), S 2 (green line), S 3 (blue line) and 
D = 9 X 10“^ (purple line). 



Figure 8 : (a) Trajectory of the center of mass (green line) and the 
shape change (red lines) of a migrating (from the left to the right) cell 
with interval St = 25 in the presence of the nonlinearity. Note that 
the linear size of this figure is about twice as large as that of Fig. EI a). 
(b)Time-evolution of A = 9/{27r) (red line), B = 02/(27r) (green 
line), and C = 9^/(27r) (blue line) for the nonlinear model. It is 
evident that the velocity direction 9 exhibits a correlated time-change 
with the angles of the deformations 92 and ^ 3 . (c) Time-evolution of 
V (red line), S 2 (green line), S 3 (blue line) and D = 9 x 10“^ (purple 
line). 
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moments is about 2 — 3 minutes. The migration (end-to- 
end) distance is about 40 /im during 15 minutes. Since the 
diameter of a cell is about 15 /rm, this means that a cell 
migrates roughly the distance of one cell size in two cycles 
of deformations. On the other hand, the average period 
of the deformation force in our model has been set about 
15 in the unit of K 2 = 1. Figure [TKa) displays the tra¬ 
jectory for the duration time 175 (about 11 cycles) during 
which the cell migrates at the distance about 10 (end-to- 
end distance). Since the radius is set as i?o = 1, that is, 
the diameter is equal to 2 , this is qualitatively consistent 
with the experiments. 


5. Discussion 


We have proposed a model system for cell crawling 
which is a set of nonlinear equations for deformations cou¬ 
pled through the migration velocity. We have introduced 
random noises for the time-evolution equations of defor¬ 
mations but not for the equation of the migration velocity 
since our basic assumption is that the internal forces cause 
crawling. When the deformation forces are constant, we 
have found that the model has a limit cycle solution. This 
oscillatory behavior is compared qualitatively with that in 
the phase field model coupled with the polarization field 
for keratocyte cells [ 2 ^ . 

In the case of time-dependent forces, they are generated 
by the coherence resonance. As a result, the deformations 
S 2 and S 3 exhibit a steep change in time. Because of the 
noise term Cz in eq. (EH), the velocity v and the phases 
62 and 63 as well as S 2 and S 3 are stochastic. However, 
this stochasticity is not strong enough to make a random 
motion like in Fig. [7Ka). Therefore we have introduced 
another set of random noises and r/n as eqs. (EH - EH) 
in the time-evolution equations. This is a technical point 
in the present theory. We do not claim that two different 
kinds of noises are involved in the motility of real living 
cells. 

For milder time-dependence, one may introduce oscilla¬ 
tory (non-stochastic) deformation forces such as 


9'^>(() = 9?>+9r (i±r^ 




( 68 ) 


(69) 


where gi^^ and gj^^ are positive constant. We 

have carried out numerical simulations of eqs. E3, EH), 
EH, EH, EH, and EH with (EH and EH with either 

$ = 7 r /4 or $ = 7 r/ 2 , and 0 < w < 27r. The results 
both for the linear case and the nonlinear case are not 
much different qualitatively from those in the excitable 
deformation forces shown in the preceding section. By 
utilizing the oscillatory deformation forces in the present 
model, the motility of human hematopoietic stem cells is 
now under investigation [37 1. 


We have emphasized that the phase difference of the de¬ 
formation forces and the relative angle of the dynamical 
variables are crucial to compare our model with experi¬ 
ments. For example, we have chosen = tt through¬ 
out the present paper. Anti-phase oscillation of S 2 and S 3 
for the constant deformation forces seems to be consistent 
with the motility of keratocyte obtained theoretically 
as shown in section 3. In section 4, the time-dependence 
of the deformation forces are chosen such that the value of 
moves counter-clockwise as shown in Fig. El or 
the phase $ in eq. (IHIHl is 0 < $ < TT. 

We have chosen the direction of deformation forces as 
eqs. (EH-(EH- It is mentioned here briefly that what hap¬ 
pens if one uses the relations (EH"(EH- I^^ this case, equa¬ 
tions for the amplitudes and the phases of deformations 
are given by 


dS2 

dt 

dt 

dS3 

dt 

dt 


-K 2 S 2 + ff^^Ht)cos(66(23 + 2T„), (70) 

-^V^sin(6d23 + 2T,), (71) 

-K3S3 + 5^^Ht) 003(6023 + 3^'„), (72) 

-^^i^sin(6023 + 3^„). (73) 


For the sake of clarity, we have written down the equations 
for the linear case without noises. If one sets = tt, the 
stationary solution of the amplitude S 3 turns out to be 
negative since K 3 and g^^\t) are assumed to be positive. 
If 'ki, = 0, i.e., 7 > 0, this difficulty is removed, but in this 
case, the cell takes the same shape as in Fig. El but mi¬ 
grates from the left to the right with the S 2 deformation fol¬ 
lowed by the S 3 deformation. Such a motion, however, has 
not been observed experimentally in Dictyostelium cells. 
See, e.g., ref. [l^ . 

We have shown that the correlation between the defor¬ 
mations and the migration velocity plays an important role 
in the cell crawling. First of all, it is noted that the re¬ 
lation EH causes a correlation between them even when 
the nonlinear coupling is absent. Secondly, when the non¬ 
linearity is present, the motion becomes more persistent. 
That is, the cell tends to migrate without large changes of 
the direction although the intensity of the random noises 
is the same as in the linear case. We note that Maeda 
et al have observed a similar dynamics experimentally in 
Dictyostelium cells Q . They have made a systematic com¬ 
parison of the motility in a starved condition and a veg¬ 
etative condition. A cell in a starved condition migrates 
at a larger velocity and elongates in the direction parallel 
to the migration direction. Its trajectory is much more 
persistent than in a vegetative condition. In the present 
analysis, we have changed the coefficients of the nonlinear 
terms to compare the linear and nonlinear cases explicitly. 
It is worth mentioning, however, that the nonlinear effects 
are enhanced for larger migration velocity by increasing 
I 7 I in eq. (15^ keeping other parameters fixed. It is an in- 
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teresting future problem to elucidate further both theoret¬ 
ically and experimentally the difference of the dynamical 
behavior between the starved and vegetative conditions. 

In the present study, we have considered only the second 
and third mode deformations assuming that higher modes 
relax sufficiently rapidly and the forces acting on them are 
weaker. Preliminary numerical investigations have been 
carried out by including the fourth mode with the relax¬ 
ation constant comparable with those of the second and 
third modes. In this case, as mentioned at the end of 
section 3, circular motions appear for time-independent 
forces in the absence of noises. On the other hand, when 
the forces of the fourth mode are time-dependent as eqs. 
(|68l) and (IS^ . spinning motion, quasi-periodic motion and 
chaotic motion have been found. We do not go into detail 
here but will publish these results separately elsewhere. 

Obviously, a remaining problem is to derive the time- 
evolution equations o, and ([3]) starting either from 
the phase field model of cell crawling, which involves the 
traction force explicitly or from the model of ac¬ 
tive droplets subjected to active stress [ 2 ^ by means of 
the reduction theory as was done for traveling domains 
in reaction-diffusion systems and for hydrodynamic 
droplets self-propelled by the Marangoni effect [s^ . How¬ 
ever, this is left for a future study. 
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